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Abstract 

We report a measurement of the Bg meson lifetime from Bg DfTX decays, where 
D~ mesons are reconstructed in the D~ — > (/)7t~ and D~ — > K*''K~ decay channels. 
From approximately 3.7 million hadronic decays recorded by the OPAL detector 
at LEP a sample is selected containing 911 ± 83 candidates, of which 519 it 136 are 
estimated to be from Bg meson decays. Fitting the distribution of the distance from 
the beam spot to the decay vertex of the D~ candidates with an unbinned likelihood 
function we measure 

^l^sJ — -L- '^-0.19-0.17 PS' 

where the errors are statistical and systematic, respectively. 
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1 Introduction 



The lifetimes of b-flavored hadrons are related to both the strength of the b quark coupling 
to c and u quarks, and to the dynamics of b hadron decay. The spectator model assumes 
that the light quarks in b and c hadrons do not affect the decay of the heavy quark, and 
thus predicts the lifetimes of all b-hadrons to be equal. For charm hadrons this prediction 
is inaccurate; the measured D"*" lifetime is approximately 2.5 times that of the D'^ and more 
than twice that of the D~ More sophisticated models predict that the differences among 
b-hadron lifetimes should be much smaller than those in the charm system, because of the 
larger mass of the b quark p|, |, §]. Bigi et al. predict a difference in lifetime between 
the B"*" and B° meson of several percent, and between the Bg and B° mesons of about 
1%. Although some assumptions in reference have been questioned by Neubert and 
Sachrajda there is agreement that the models predict a difference between the B^ and 
B° lifetimes of the order of ±1%. 

The first measurements of the Bg lifetime ^ were made using correlated D~-lepton pairs 
that primarily result from the semileptonic decay of the B^. However, the small Bg semilep- 
tonic branching ratio limits the statistical power of this channel. More inclusive techniques 
can be used to circumvent this limitation. 

We present a new measurement of the lifetime of the B^ meson in which only the D~ mesons 
are reconstructed. The two decay channels used are:Q 

B° ^ D- X B° ^ D- X 

U K*oK- and U ^vr" 

U K+TT- U K+K-. 

The Bg meson is not the only source of Dg" mesons. Significant numbers of D~ mesons are 
produced in the decays of other b-hadrons as well as in 1? — > cc events. The level and shape 
of this background is evaluated using Monte Carlo data and measured branching fractions. 
The background from random tracks combination is evaluated using the side band technique. 

To measure the Bg lifetime, we reconstruct the decay vertex of the D~ and determine the 
distance from the beam spot to this point. For D~ mesons that result from the decay of 
b-hadrons, this "decay length" has contributions both from the b-hadron and Dg* decay 
lengths. The bias arising from the angle between the Bg meson and the D~ meson direction 
is very small and is not affecting the result of this analysis. An unbinned maximum likelihood 
fit is performed using the reconstructed decay lengths, their errors and the D~ momentum 
to extract the mean Bg lifetime. The following sections describe the OPAL detector, the 
selection of D~ candidates, the vertex topology of the events, the determination of the Bg 
decay length, the estimation of the Bg energy, the lifetime fit, the results, and the systematic 
errors. 



^Charge conjugate modes are always implied. Also, unless otherwise noted, K and tt always refer to 
charged particles. 
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2 The OPAL Detector 



The OPAL detector is described in reference 0. The central tracking system is composed of a 
sihcon microvertex detector, a precision vertex drift chamber and a large volume jet chamber 
surrounded by a set of chambers to measure the 2;-coordinate of tracks (2;-chambers)0. These 
detectors are located inside a solenoid. The detectors outside the solenoid consist of a time- 
of-flight scintillator array and a lead glass electromagnetic calorimeter with a presampler, 
followed by a hadron calorimeter consisting of the instrumented return yoke of the magnet, 
and several layers of muon chambers. Charged particle types are identified by their specific 
energy loss, dE/dx, in the jet chamber. Further information on the performance of the 
tracking and dE/dx measurements can be found in reference |p. 



3 Dg Candidate Selection 

This analysis uses data collected during the 1991-1995 LEP running periods at center-of- 
mass energies within ±3 GeV of the mass. After the standard hadronic event selection 
and detector performance requirements, a sample of 3.7 million events is selected. Charged 
tracks and electromagnetic clusters not associated with a charged track are grouped into 
jets using the JADE EO recombination scheme |]10[ with a ?/cut value of 0.04. Tracks from 



identified secondary vertices, A and Kg decays or 7 conversions, are excluded from the 
candidate selection. 



Simulated event samples were generated using the JETSET 7.4 Monte Carlo program [11 



together with a program to simulate the response of the OPAL detector H]. The Monte 



Carlo sample includes approximately 4 million simulated multihadronic Z° decays and one 
million Z° bb decays (the equivalent of about 4.5 million multihadronic decays). In ad- 
dition, three special Monte Carlo samples were generated in which each event contains at 
least one D~ decaying in the channels of interest. The parameter optimisation used in this 
simulation is described in [rT|. For each of the following decay channels, 20000 events were 



generated: — > X, Dg from b-hadrons other than D^ decays and Dg from Z — >■ cc. 



3.1 D3- K+K TT" selection 

The D~ meson is reconstructed in the decay chains D~ — > K*°K~ in which the K*'^ decays 
into a K+TT", and D~ 07r~ where the subsequently decays into K^K~. 

Tracks forming the D~ candidates are required to be in the same jet and to have the appro- 
priate charge combination. At least two of the three candidate tracks are required to have 
good 6 measurements either from the z-chambers or from a measurement of the track end- 
point as it exits the main jet chamber. Similarly, to reject poorly reconstructed candidates, 
at least two of the three tracks are required to have hits in the silicon microvertex detector. 

^The coordinate system is defined such that the z-axis follows the electron beam direction and the x- 
y plane is perpendicular to it with the a;-axis lying horizontally. The polar angle 9 is defined relative to the 
+z-axis, and the azimuthal angle (f> is defined relative to the +j:-axis. 
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To reduce the combinatorial background, the tracks forming the D~ are subject to particle 
identification requirements. For candidate pion tracks, the probability for the measured 
dE/dx value to be consistent with the pion hypothesis is required to be greater than 1%. 
For candidate kaon tracks, if the observed energy loss of a kaon candidate is less than the 
mean dE/dx expected for a kaon, the probability of consistency with the kaon hypothesis 
is required to be greater than 1%, and greater than 3% otherwise. If both kaon candidates 
have energy losses greater than the mean dE/dx expected for a kaon, the product of the two 
dE/dx probabilities is required to be greater than 0.02. These tighter requirements reduce 
the background from pion tracks, for which the mean dE/dx value is above that of kaons. 
Furthermore, one of the kaon candidates must satisfy a pion-rejection criterion by having an 
observed dE/dx less than that expected for a pion, and a probability of less than 10% that 
the dE/dx is consistent with a pion hypothesis. 

For the K*°K~ mode, the invariant mass of the Ktt combination is required to satisfy 0.865 < 
rriKn < 0.925 GeV. In the (pn mode, the width of the K+K" peak is dominated by detector 
resolution and the K"'"K~ invariant mass is required to satisfy 1.010 < mxK < 1.030 GeV. 
The momenta of the kaons are required to be greater than 2 GeV and that of the pion must 
exceed 1 GeV. 

To reduce further the D~ K*^7r~ background in the K*°K~ mode, the kaon candidate 
originating directly from the Dj7 decay is subject to a tighter, 5%, dE/dx requirement. 
Also, in this mode both kaons must meet the pion-rejection criterion described above. 

We require that the D~ momentum divided by the beam energy, Xj-,-, is greater than 0.20 to 
reduce random track combinations, and less than 0.60 to reduce the D~ contribution from 
Z° — > cc events. 

The differences between the angular distributions of D~ decays and those of combinatorial 
background are exploited to enhance the signal purity as follow. The D~ is a spin-0 meson 
and the final states of both decay modes consist of a spin-1 (0 or K*'') meson and a spin-0 
(vr" or K~) meson. The D~ signal is expected to have no dependence on cos6p, where 6p 
is the angle in the rest frame of the D~ between the spin-0 meson direction and the D~ 
direction in the lab frame. However, the cos 6p distribution of random combinations peaks 
in the forward and backward directions. It is therefore required that | cos6'p| < 0.8 for both 
modes. The distribution of cos 6^, the angle in the rest frame of the spin-1 meson between the 
direction of the final state kaon from the decay of the spin-1 meson and the D~ direction, 
is proportional to cos^ for D~ decays. The cos^„ distribution of the random K+K~7r~ 
combinations in the data is, however, approximately uniform. Therefore it is required that 
I cos^^l > 0.6 (0.4) for the K*K (07r) mode. 



3.2 Decay length determination 

The D~ decay vertex is reconstructed in the x-y plane by fitting the D~ candidate tracks to 
a vertex. To reject very poorly reconstructed vertices, the of the vertex fit is required to 
be less than 10 (for 1 degree of freedom). 

The beam spot position is measured using charged tracks in the opal data with a technique 
that follows any significant shifts in the position during a LEP fill |l^. The intrinsic width 
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of the beam spot in the y direction is about 8 /im. the width in the x direction is directly 
measured using events and varies between 100 /im and 160 /zm, depending on the LEP 

optics. 

The distance from the D~ decay vertex to the beam spot is determined in the x-y plane. This 
distance is converted into three dimensions using the polar angle of the K"*"K^7r~ momentum 
vector. Typical decay length errors are about 300 yum, with only a small contribution coming 
from uncertainties in the position of the interaction point within the beam spot. Rejecting 
candidates with decay length errors greater than 3 mm helps to reduce the effects of poorly 
measured tracks. 



3.3 Dg selection results 



The K"'"K~7r~ invariant mass distribution for all candidates that pass the selection is shown 
in figure |1|. A likelihood fit to the measured K+K^tt^ invariant mass of the D~ candidates is 
performed separately for the two decay channels. The K+K"??" mass distribution is param- 
eterized as the sum of a quadratic term to account for random combinatorial background 
(which also is observed to describe accurately the background in the simulated data sam- 
ples), a Gaussian function which describes the mass peak of the reconstructed D~ signal, 
and another Gaussian to account for a K+K^vr^ contribution. This last term has a 

mean fixed to the D~ mass of 1869.3 MeV ll| and the width constrained to be the same as 
that of the D~ signal. This second Gaussian is included to avoid biasing the estimate of the 
combinatorial background. The results of this fit to the K"'"K~7r~ invariant mass spectra are 
shown in table A total of 911 ± 83 D~ candidates are found. 



source 


07r 


K*K 


number of Dg" candidates 


629 ± 64 


282 ± 53 


fitted D~ mass 


1.966 ± 0.002 GeV 


1.959 ±0.004 GeV 


fitted D~ sigma 


18 ± 2 MeV 


21 ± 3 MeV 


background fraction 


0.75 


0.86 


total over 80 bins 


90.0 


80.0 



Table 1: K+K^vr^ invariant mass spectra fitting results. The fitted mass is in agree- 
ment with the value of 1.9685 it O.OOOGMeV from The "background frac- 
tion" is defined as the fitted number of combinatorial background candidates 
divided by the total number of candidates within 2a of the fitted D~ mass. 



3.4 Composition of the signal 

The Bg meson is not the only source of D~ mesons. Significant numbers of D~ mesons are 



s 

7Q 



produced in other b-hadron decays as well as in Z cc events, collectively referred to as 
'physics background'. In what follows, the contributions of these three components to the 
observed D~ spectra are estimated. 

The fractions of D~ mesons produced from each of these sources are extracted using the 



results in reference 0]. The measured production rates separated into flavour are 
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Figure 1: Results of the standard selection. Top: K+K~7r~ invariant mass distribution 
for (pTT combinations along with the fitted distribution. Bottom: K"'"K~7r~ 
invariant mass distribution for K*K candidates, along with the fitted distri- 
bution. The dotted lines indicate the region within which candidates are used 
in the lifetime fit. The arrow indicates the D~ peak position. 
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/b = Tbb/rhad ■ /(b ^ D-) • Bi{D- ^ 07r-) = (0.166 ± 0.018 ± 0.016)% and = T,-JT,,^^ ■ 
/(c D-) • Br(D- ^ 07r-) = (0.056 ± 0.015 ± 0.007)%. Thus, /b/(/c + /b) = (75 ± 9)% 
of produced D~ mesons are from Z'^ — bb events and the remaining (25 ±9)% are from 
Z" — »• cc events. 

The fraction of D~ mesons from b-hadrons in our signal is estimated using the following 
production rates: /(b ^ 6°) = 0.112;[!:[]J^ and /(b ^ B) = 0.378 ± 0.022 |I[ (where 'B' 
is either B"*" or B°) . Assuming that b-baryons decay to Dj7 mesons with the same branching 
fractions as the non-strange b-mesons (0.086 ± 0.016 |]1|) and using the inclusive branching 
ratio of Bg to D^, measured to be 0.87 ±0.31 [jl|, we estimate that (56 ±11)% of D~ mesons 
from b-hadron decay are from B^ . Thus, (42 ± 10)% of the D~ mesons come from B^ decays, 
(33 ±9)% come from other b-hadron decays and (25 ±9)% from Z° — cc events. Assuming 
the production rates of B"*" and B'' are the same and that they have equal branching ratios 
to D~, each of these non-strange b-mesons accounts for (14 ±4)% of the D~ production with 
the remaining (5 ± 2)% of D~ produced from b-baryons. 

The special Monte Carlo samples described above were used to determine the contribution of 
each of these channels i the reconstructed sample. The ratio of the efficiency to reconstruct 
a Dg" meson from b-hadron decay other than B^, divided by the efficiency to reconstruct a 
Dj7 meson from a B^ decay is 0.82 ± 0.02, where the error is due to the limited statistics of 
the simulated data samples. The principal reason that this ratio is less than unity is that 
the momentum spectrum of D~ mesons from the decay of b-hadrons other than B^ is softer 
than that of D~ mesons from B^ decay. The ratio of the efficiency to reconstruct a D~ 
meson from a Z° — cc event divided by the efficiency to reconstruct a D~ from a B^ decay is 
0.67 ± 0.01. The upper cut on the scaled energy of the D~ (xj-,- < 0.6) preferentially rejects 
D~ from Z° — cc, which tend to have higher momentum, the sources of D~ production are 
summarized in table 0. 



source 


fraction of produced Dg 


fraction of reconstructed Dg 


Z" cc 


25 ±9 


17 ±6 


B+ 


14 ±4 


11 ±3 


BO 


14 ±4 


11 ±3 


b-baryons 


5±2 


4±2 




42± 10 


57± 14 



Table 2: Estimated D~ signal composition. The errors include those from the measured 
branching ratios and the statistical uncertainty from the Monte Carlo samples 
used to estimate the relative reconstruction efficiencies. 

Monte Carlo events were used to study the background from events where the three K+K~7r~ 
candidate tracks come either from the same fully reconstructed charm hadron for which a 
pion or a proton has been mis-identified as a kaon, or from a partially reconstructed charm 
hadron. The resulting K+K~7r~ invariant mass distribution in the region around the D~ mass 
is similar to that of the combinatorial background. Such events are therefore, considered to 
contribute to the combinatorial background. 

Since (57± 14)% of the reconstructed D~ meson decays are calculated to be from Bg decays, 
519 ± 136 D~ candidates are attributed to Bg decay from the 911 ± 83 candidates resulting 
from the fit to the K+K~7r~ invariant spectra. 
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4 The Lifetime Fit 



To extract the lifetime from the measured decay lengths, an unbinned maximum likelihood 
fit is performed using a function that accounts for both the D~ signal and the background 
components of the sample. In the part of the likelihood function describing events from 
b-hadron decays, the observed decay lengths depend on the Bg lifetime. 



The form of the likelihood function for the candidates in the 
described in terms of the probability for observing a combined D~ 
given a measurement error a}^, the momenta of the D~ ""'^ 
these mesons. 



signal from B^ decays is 
and Bg decay length, L\ 
and Bg, and the mean lifetimes of 



The likelihood function has components which describe the different sources of D~ mesons 
in the signal and in the combinatorial background. This follows closely the method used in 
previous opal B^ lifetime analyses 0. 

The likelihood function which accounts for D~ mesons from B^ decays is constructed as a 
convolution of two exponential functions to describe the D~ and B^ decay lengths, convoluted 
with a function to describe the probability of having a particular Bg momentum (pbs) and 
Gaussian functions to describe the measured decay length resolution. This can be expressed 
as: 



rO 



dl / dpB,7^(L' I /, al, si, S2, /2) B{pb, \ p^J V{1 \ tbo, rD-,PB,,PD, 

JO s s s s 

where p-^X'^'^^) = 45GeV and IZ is given by 

n{U I /, ai, Si, S2, /2) = (1 - /2) <3{L^ I /, ai, si) + h G{L' I /, al, s^) . 



(2) 



The function } is a Gaussian function which describes the probability to observe a decay 
length, L*, given a true decay length I and the measurement uncertainty cr}^ and scale factors 
on this error, si and S2- Two scale factors are employed to describe both the majority of 
tracks for which the measured decay length uncertainty is a good estimate and the small 
fraction, /2, of mis-measured candidates in which the measured decay length uncertainty 
is an underestimate. The scale factors Si and S2, as well as the fraction of mis-measured 
candidates, /2, are free parameters in the lifetime fit. B is the probability of a particular Bg 
momentum for an observed D~ momentum,pDs- This probability is determined from Monte 
Carlo events by forming distributions of the ratio PbJpDs- Six such distributions are formed, 
depending on the value of p^^, since at higher values of p^^ the range of potential values of 
PBs is more tightly constrained than for lower pd^ candidates. V is the probability for the 
Dj7 to decay at a distance / from the e~^e~ interaction point, given Bg and D~ lifetimes t^o 
and Tj-)- and momenta ps^ and pl^^. This function is constructed as: 





-I' ■ rriB, 




/ exp 
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where tboPBs/''^Bs and Tj^-pdJitld^ are the mean decay lengths for the given momenta (pbs 
and pdJ mean hfetimes (tbo and Tj-,-) and masses {m-B^ and ^dJ- 

Similar functions are employed for the other D~ meson sources. For D~ mesons from other 
b-hadron decays (B+, B° and b-baryons), the world average lifetime for each species of 
b-hadron [|l| is used in the likelihood function and a slightly different boost function, B, is 
employed. For D~ mesons produced in Z° cc events, the function P is a single exponential 
function with the decay constant Tj-,-pds/'"^Ds- The total likelihood function containing the 
contributions for all sources of D~ mesons, is formed by combining the likelihood functions 
for each of the sources of D~ mesons with the fixed fractions listed in table ^.0 This likelihood 
is written as: 

Cf^{V\T^o,alpi,J = (1 - - /b+ - /bo - /Aj£f + 

fc-c ■ Cf + /b+ ■ Cf + /bo ■ Cf + ■ . (4) 

The functional form of the likelihood function, used empirically to parameterize the 

combinatorial background, is composed of an exponential convolved with the boost func- 
tion, B, a fraction with no lifetime and the same double-Gaussian resolution function as the 
signal. This is expressed as: 

/:--b(LMr,g,/o,ai,pb.,Si,S2,/2) = 

dl dpBML' I I, al si, S2, /2)S(PB. I PnjrtgiL \ n^, f,p^J , (5) 

where the parameters describing the resolution, si, S2 and /2, are the same as used in the 
likelihood terms that describe the D~ signal and 

+ fS{l) . (6) 

The fraction of background with no lifetime, /°, as well as the characteristic lifetime of the 
background, Tbg, are free parameters in the fit. 

The combinatorial background in the event sample is taken into account by fitting for it 
simultaneously with the D~ signal. The probability that a candidate has arisen from a 
combination of background tracks, jF™™^(mj), is determined as a function of the K"'"K~7r~ 
invariant mass of each candidate, rrii, using the results of the fit to the K"'"K~7r~ invariant 
mass spectrum. 

Thus, the full likelihood for candidate i is: 

Ci{V I TBO, al,^,, = [1 - ^--''(m,)] ■ Cf^ + ^-™^(m,) ■ CT""' . (7) 

In total six parameters are free in the fit: the B^ lifetime, the parameters describing possible 
scale factors on the decay length error (si, $2 and /2) and the parameters describing the 
combinatorial background (/° and Thg). 

■^The momentum dependence on The fractions /^g, /b+ , f^o and has been neglected. This omission 
was studied in simulated data and found not to produce a noticeable bias in the resulting lifetime. 



Vbg{l I Ug, f,PBj = (1 - / ) exp 
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The lifetime fit uses the 10633 events found in the region from 50 MeV below the fitted D~ 
mass to 200 MeV above it (see figure |I]). From studies on simulated data, it is found that 
the lifetime properties of the combinatorial background in this region accurately reflect the 
combinatorial background around the D~ mass. Furthermore, this avoids the region below 
the D~ mass in which the number of K'^K^n^ candidates from D~ decays and contribu- 
tions from other D meson decays (e.g. reflections and partially reconstructed decays) are 
potentially significant. This fit finds T(Bg) = 1.72lQ^gps, where the error is statistical only, 
the values of the other free parameters of the fit are in table |[ The fitted values of the 



parameter 


fit results 


/2 


0.07 ±0.01 


Si 


0.88 ±0.01 


S2 


4.69 ±0.26 


f 


0.77 ±0.01 




1.61 ±0.07 



Table 3: Final values of the free parameters in the fit. 

parameter are describing the decay length resolution are consistent with the understanding 
of the OPAL tracking performances. The decay length distributions are shown in figure ^ 
separately for candidates with K"'"K~7r~ invariant mass within 50 MeV of the fitted D~ mass 
(the "signal region"), and for those candidates outside this mass window (the "sideband 
region"). These illustrate the quality of the fit in regions dominated by D~ signal events and 
by combinatorial background, respectively. The curves in figure ^ represent the sum of the 
decay length probability distributions for each event. Using the 42 bins that are expected to 
contain at least five candidates (as predicted by the lifetime fit), a total of 54.7 is found 
for the sum of the signal and sideband decay length distributions. For the positive decay 
length bins, a total of 24.0 for 29 bins is observed. These values and plots shown in 
figure indicate that the fitted functional forms provide a good description of the data for 
both signal and background. It should be stressed that the fit is to the unbinned data. 



5 Checks of the Method 

A number of different checks have been made to investigate potential biases in the method 
of selecting and fitting the signal. 

5.1 Potential bias in the selection and fitting procedure 

Tests were performed on several samples of simulated data to check for biases in the selection 
and fitting procedures. The first tests involved a toy Monte Carlo program which generated 
decay length data for the signal D~ decays and combinatorial background. For each D~ 
signal candidate from a decay, this simulation generated and D~ decay times from 
exponential distributions with the means set to known values. The Bg and D~ momenta 
were chosen from a spectrum based on the full Monte Carlo simulation. The B^ and D~ 
decay lengths were then calculated and combined to give the true candidate decay length. 
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Decay Length (cm) 



Figure 2: Top: decay length distribution within a D~ mass region of zb50MeV around 
the fitted D~ mass (signal region). The single-hatched area represents the 
contribution from physics background, the unhatched area the combinatorial 
background and the cross-hatched area is the contribution from decays. 
Bottom: decay length distribution for events outside the mass region (side- 
band region), namely from 50 to 200 MeV above the D~ mass. 
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which was then smeared by a resolution function, physics and combinatorial backgrounds 
were generated through a similar procedure. Many fits were conducted over wide ranges of 
Bg lifetimes with different levels and parameterizations of the backgrounds. The result of 
these studies shows that any bias in the fitted 3° lifetime is less than 0.5% and that the 
statistical precision of the fit to data is consistent with the sample size and composition. 

To verify that the D~ selection does not bias the reconstructed sample, a lifetime mea- 
surement was made from 20000 — * DsT^ Monte Carlo decays into the two channels of 
this analysis, using a lifetime of 1.60 ps. The mean lifetime of the selected D~ sample 
was 1.64 ± 0.04 ps, consistent with the expectation that there is no bias from the selection 
procedure. The lifetime obtained by fitting this same sample was 1.65 ± 0.05 ps. 

To investigate the effects of the combinatorial background on the lifetime fit, the same 
selection and fitting procedure has been applied on a Monte Carlo sample of 4 million 
multihadronic zq decays. Due to the choice of branching ratios used to produce this sample 
of simulated events, there are fewer reconstructed D~ signal candidates than we observe in 
the data. The fitted lifetime has been found to be 1.73 ± 0.29 ps, which is consistent with 
the generated B^ lifetime of 1.6 ps within the statistical power of this sample. If the signal 
events from the 20000 B^ — > D^Tx decays described above are added to the simulated data, 
the resulting sample is of similar purity to that found by the tight neural net selection. The 
fitted lifetime of this pure sample is 1.65 ± 0.07 ps, again in good agreement with the true 
lifetime of the sample. 

The lifetime fit has also been repeated as for the standard result, except that the (pn and K*K 
channels are fitted separately. The results are 1.53 ± 0.23 ps and 2.14±0.40ps, respectively, 
consistent at the level of 1.3 standard deviations. 

5.2 Use of the tight (neural network) selection 

As a check, a much tighter selection was developed which employs a neural network to 
reject significantly more combinatorial background, thereby producing a much purer D~ 
signal. However, this also results in a rather significant loss of signal, and as such is not 
as statistically powerful for the B^ lifetime determination. This artificial neural network 
uses 16 kinematic and particle-identification quantities, including b-tagging information in 
the hemisphere opposite to the D~ candidate. The K~^K~7t~ invariant mass distributions 
for all candidates which pass this selection are shown in figure |^. The likelihood fit to the 
measured K"'"K~7r~ invariant mass of the D~ candidates is performed as above, resulting in 
a smaller signal of 232 ± 29 (104 ± 17) candidates with a much reduced background fraction 
of 0.54(0.61) for the 07r (K*K) mode. 

The sample of D~ candidates found by the tighter neural network selection has also been 
used to fit for the B^ lifetime. This sample gives a value of 1.69 ± 0.27 ps for the B^ lifetime 
and is consistent with the more precise result from the analysis of section 4. The differ- 
ence in statistical precision between the two fits is in agreement with expectation from toy 
Monte Carlo studies, given the relative size of the D~ signal and the level of combinatorial 
background in each sample. 

The results of all of the tests show no indication of a significant bias in the selection and 
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1.6 1.8 2 2.2 2.4 

KKtt invariant mass (GeV) 

Figure 3: Results of the tight (neural network) selection, used as a cross-check of the 
analysis. Top: K+K~7r~ invariant mass distribution for ^tt combinations along 
with the fitted distribution. Bottom: K~^K~tt~ invariant mass distribution for 
K*K candidates, along with the fitted distribution. The dotted lines indicate 
the region within which candidates are used in the lifetime fit for the check. 
The arrow indicates the d~ peak position 
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fitting procedures. 



6 Evaluation of Systematic Errors 

Systematic errors arise from tlie level, parameterization and source of the background, the 
potential bias from the selection and fitting procedure, the boost estimation method, the 
beam spot determination and possible tracking errors. The systematic errors are summarized 
in table |. 

Combinatorial background 

We consider the effects of both the level of the combinatorial background, as determined by a 
fit to the K+K^7r^ mass spectrum, and the effective lifetime of this background, as estimated 
from the candidates in the sideband region. The systematic error due to the level of this 
background is determined by repeating the calculation of the event-by-event probability that 
an event is combinatorial background by changing the estimated size of the D~ signal by the 
statistical uncertainty from the fit to the invariant mass spectrum. This yields a variation in 
the Bg lifetime of ±0.08 ps. The width of the K^'^K^tt^ mass region from which candidates 
are selected for use in the lifetime fit has also been varied. This was done by selecting 
candidates in regions extending from 150 to 350 MeV above, and from 25 to 75 MeV below, 
the D~ mass. Another check included using a sideband from 200 to 350 MeV above the D~ 
mass in place of the standard sideband from 50 to 200 MeV. These change the fitted 
lifetime by ±0.07 ps which is assigned as a systematic error. Tests were also conducted using 
the toy Monte Carlo which indicated that the level of these observed changes is consistent 
with the expected 0.05 ps uncertainty due to statistical fluctuations in the sideband sample; 
however we conservatively retain the observed variation as a systematic uncertainty. 

Several alternative parameterizations describing of the decay length distributions of the com- 
binatorial background have been investigated. For example, we have included an exponential 
on the negative decay length side, in place of the second Gaussian function, to describe those 
events that are significantly mis-measured. The resulting Bg lifetime is ±0.02 ps higher than 
the standard result. In another check the second wider Gaussian was used only for the com- 
binatorial background term in the likelihood function, changing the lifetime by ±0.05 ps. 
Leaving out this second Gaussian altogether decreased the lifetime by 0.10 ps, although the 
quality of this fit is significantly worse, as illustrated by an increase in of about 300 for 
40 bins, c.f. 54.7 for 42 bins with the default parameterization. Consequently this last case 
is not considered as a systematic error. In no case do these alternative parameterizations 
significantly improve the quality of the B^ lifetime fit, and an error of ±0.05 ps was assigned 
to cover such effects. 

We assign a total systematic error due to the combinatorial background, parameterization 
and source, of ±0.12 ps. 

Physics background: sources of D^T mesons 

The physics background composition has already been discussed. Varying the fraction of D~ 
mesons from Z° cc events over the range given in Table produces a change of Il^'H ps 
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in the lifetime. The uncertainty due to the D~ fraction from b-hadrons other than 
is evaluated in the same way. The observed change on the B^ lifetime is lo;o3 ps. These 
variations also change the statistical error on the fitted B^ lifetime by up to ±0.03 ps, and 
this systematic effect is not considered further. 

The Bg lifetime dependence on the uncertainty of the b-hadron lifetimes has been measured 
by varying the b-hadron lifetimes within the errors quoted in reference [0] , assuming, conser- 
vatively, that the individual lifetimes are fully correlated with each other. The B^ lifetime 
changed by ±0.02 ps, which is included as a systematic error. The D~ lifetime has also been 
varied in the fit, within the errors quoted in reference The effect on the measured B^ 
lifetime is ±0.01 ps. This variation is also taken as a contribution to the systematic error. 

The D~ momentum spectrum also depends on the fragmentation in TP cc and Z° bb 
events. Changes in the fragmentation affect the composition of the D~ signal through changes 
in the efficiencies for the charm and bottom contributions to the D~ signal (the effect of the 
momentum spectrum of the D~ in TP bb events on the estimation of the B^ momentum 
in the lifetime fit is discussed below). In Z° — >■ bb events, we have varied the average b 
hadron energy by the measurement errors [T^ to yield a variation in the observed lifetime of 



±0.01 ps. Similarly, the momenta of the D~ mesons produced in Z° cc events was varied 
according to the average energy measured for non-strange mesons in Z" cc events (on 
the assumption that there is little difference in fragmentation amongst the various charm 
mesons) fl^, producing a variation of ±0.01 ps in T(Bg) 



In the case of D~ mesons coming from the decay of b-hadrons other than B^, the D~ may be 
produced in either a two-body mode (e.g. B Dg^D) or a multi-body final state where one 
or more light particles are produced. The two-body decay fraction of B — >^ D~X has been 
measured to be 0.56 ± 0.10 |]16|. In determining the relative efficiency of these D~ mesons 



from b-hadrons other than B^, we have already corrected our simulation to this two-body 
fraction. Assuming b-baryons decaying to D^T mesons in two- or multi-body states have 
the same fractions as B+ and B° mesons, the efficiency of selecting D~ mesons from these 
decays has been evaluated by varying the two-body fraction over the range 0.46 to 0.66 and 
re-evaluating the relative efficiency of this source of D~ with respect to the D~ which arise 
from a B^ decay. This produces a variation of ±0.01 ps in the fitted B^ lifetime. 

Thus we assign a total error due to these other sources of D~ mesons of '^^^2 PS- 

Boost estimation 

The energy spectrum of the Monte Carlo Bg events used to estimate the momentum of the 
Bg, given the observed D~ momentum, can also affect the resulting lifetime. This effect is 
not large because the scaling used to estimate the B^ momentum from the measured D~ 
momentum is correlated with the D~ momentum itself. We have varied the average B^ 



energy by the measured errors on the average b-hadron energy to yield a variation in 
the Bg lifetime of ±0.01 ps.0 The effect of a 2.0 MeV uncertainty in the mass of the Bg [Q 
was found to produce a change of less than 0.01 ps in the Bg lifetime. 

''Note that the uncertainty due to the bottom and charm hadron energy spectra affects the Bg hfetime 
both through the boost estimation and through the D~ sample composition. When combined, these two 
contributions are added hnearly. 
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Beam position and size 



The average intersection point of tlie LEP beams in OPAL is used to estimate the production 
vertex of the 6° candidates. The sensitivity of r(Bg) to the assumed position and size of 
the beam spot has been evaluated as in reference 0, resulting in a variation in the fitted 
lifetime of no more than 0.01 ps, which has been assigned as a systematic error. 

Detector alignment 

The effects of alignment and calibration uncertainties on the 3° lifetime are estimated from a 
detailed study of 3-prong r decays |T3[ 
on t(B°). 



These uncertainties lead to an uncertainty of 0.01 ps 



source 


uncertainty (ps) 


combinatorial background 


±0.12 


physics background 


+0.13 
-0.12 


uncertainty in boost 


±0.01 


beam spot 


±0.01 


alignment errors 


±0.01 


total 


+0.18 
-0.17 



Table 4: Summary of systematic errors on the Bg lifetime. 



Combining the systematic errors in Table we find r(Bg) = 1.72;to;i9(stat)tQ;J7(syst) ps. 



7 Conclusion 



A sample of D~ decays has been reconstructed in which the D~ has decayed into K~^K~n~ 
through either the (f)TT~ or K*''K~ channels. From 3.7 million hadronic Z" events recorded 
by OPAL from 1991 to 1995, a total of 911 ± 83 such candidate decays have been found, of 
which about 57% are expected to be from decay. The lifetime is found to be 

r(B0) = 1.72t°:?°(stat)l°:J?(syst)ps, 

a result consistent with the measured B'^ lifetime and other B^ measurements [|l|. This result 
is also in agreement with current theoretical expectations ||^, ^. 
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